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Abstract 

Using Lindblad approach to study decoherence of quantum systems, we study the decoherence and 
decay of entangled states, formed by two basic states of a chain of thee qubits. We look on these states 
for a possible regular dependence on their decay as a function of their energy separation between the 
basic states under different type of environments. We found not regular or significant dependence on this 
energy separation for the type of environment considered . 
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1 Introduction 

In the real world (quantum or classical) the interaction of the system with the environment is unavoidable. 
In principle, one could study the unitary evolution of the whole system, quantum plus environment plus 
quantum-environment interaction, but this represent a many bodies problem which is unsolvable within any 
picture of the quantum mechanics. The most used approach to study this phenomenom is to use the matrix 
density approach for the whole system and to make the trace over the environment variables [2 H112 [5]. 
The resulting density matrix is called ” reduced density matrix”, and its associated non-unitary evolution 
equation is called ’’master equation.” This equation is phenomenological where a dissipative and diffusion 
parameters are defined, and they are responsible of the decay behavior of the non diagonal elements of the 
reduced density matrix. This phenomenon is called ” decoherence” because is related with the disappearance 
of the interference terms of the product of the quantum wave function mm and many think that this 
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decoherence effect is closely related with the appearance of the classical world [a 1 [7]. In most of the 
approaches, the positiveness and trace equal to one are kept as principal condition for the reduced density 
matrix. The best known mathematical approach which kept these conditions was given by Lindblad cni, 
which gave an abstract general non unitary evolution equation for the reduced density matrix. The master 
equation is different when dealing with continues systems (quantum Browning motion, for example) [TJ [ 2 ] 
or discrete quantum systems (spin system) [TT]. One of the used approaches for quantum discrete system is 
described in [T^] , and we will use this approach for our study of decoherence of entangled states built up with 
two states of three qubits in a quantum computer model of a linear chain of three paramagnetic atoms with 
nuclear spin one half m- In this work, we are interested in determine the decoherence of several entangled 
states formed by two states of three qubits, and we will use the above mentioned Markovian-Lindblad master 
type of equation [Ilia- On the other hand, even this model for solid state quantum computer has not been 
built, it has been very useful for theoretical studies about implementation of quantum gates and quantum 
algorithms m min US] which can be extrapolated to other solid state quantum computers. The main idea 
is to explore the possible sensitivity of the decay of an entangled state with respect the difference energy of 
its two states involved, we establish the four cases to be considered with the quantum-environment system: 
independent environment interaction, pure dephasing interaction, correlated dissipation interaction, and 
dephasing correlated interaction. The analytical dynamical systems of the reduced density matrix elements 
are obtained for these cases, and the results of the analytical and numerical simulations are presented. 


2 Hamiltonian of the chain of nuclear spins 


Following Lloyd’s idea consider a linear chain of nuclear spin one half, separated by some distance and 
inside a magnetic in a direction z, B(z) = (0, 0, i?o(- 2 )), and making and angle 0 with respect this linear 
chain. Choosing this angle such that cos 9 = 1/ \/3, the dipole-dipole interaction is canceled, the Larmore’s 
frequency for each spin is different, ujk = ^BQ(zk) with 7 the gyromagnetic ratio. The magnetic moment of 
the nucleus flk is related with its spin through the relation fik = 78 ^, and the interaction energy between 
the magnetic field and magnetic moments is Hint = — ■ B(zfc) = — Xa, If in addition, one has 

first and second neighbor Ising interaction, the Hamiltonian of the system is just [13] 


Hs 


N j N-1 


2J' 


N-2 


E qz qz 


( 1 ) 


where N is the number of nuclear spins in the chain (or qubits), J and J' are the coupling constant of the 
nucleus at first and second neighbor. Using the basis of the register of N-qubits, {|^^,..., ^ 1 )} with = 0,1, 
one has that 5'^|C/c) = Therefore, the Hamiltonian is diagonal on this basis, and its eigenvalues 
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are 


k^l k^l k^l 

3 Interaction with the environment 

Consider now that the environment is characterized by a Hamiltonian He and its interacting with the quantum 
system with Hamiltonian Hg. Thus, the total Hamiltonian would he H = Hg + He + Hge, where Hge is the 
part of the Hamiltonian which takes into account the interaction system-environment, and the equation one 
would need to solve, in terms of the density matrix, is [Hill 

= [H,pt], (3) 

where pt = pt{s, e)is the density matrix which depends on the system and environment coordinates. The 
evolution of the system is unitary, but it is not possible to solve this equation. Therefore, under some 
approximations and tracing over the environment coordinates EBa-, it is possible to arrive to a Lindblad 
type of equation for the reduced density matrix p(s) = tre{pt), 

= [Hg,p] + - Iv^V.P - (4) 

where Vi are called Kraus’ operators. This equation is not unitary and Markovian (without memory of the 
dynamical process). This equation can be written in the interaction picture, through the transformation 
p = UpW with U = as 

(5) 

where £(p) is the Lindblad operator 

( 6 ) 

with V = UVU^. The explicit form of Lindblad operator is determined by the type of environment to consider 
m at zero temperature. In this work we consider dissipation effects and defacing. So, the operators can be 
Vi = S~ (for dissipation), Vi = (for defacing), and ji/ih is the coupling constant with the environment. 
In this way, one considers the following cases: 

(a) Independent : In this case, each qubit of the chain acts independently with the environment, and one 
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has local decoherence of the system. The Lindblad operator is 


N 


= - ^tSuP - pstsu) 


where 5^ and 5^, are the ascend and descend operators such that 


sf = usfu^ = 


( 7 ) 


( 8 ) 


where has been defined as 

(ik =Wk + ■^(S'fe+I + S'fc-i) + -^(<S'fe+2 + Sk-2)- (9) 

(b) Correlated independent : Each qubit interact with the environment but its effect is felt by the other 
qubits, that is, the type of interaction is nonlocal with a collective interction between qubits and environment. 
The Lindblad operator is 

1 ^ 

= -fiJ2 - S^SkP - pS^Sk)- (10) 

where ^jk are the coupling constant between qubits and environment, with = ^kj and = y* . 


(c) Dephasing : There is not interchange of energy between qubits and environment, only decoherence 
is presented where the non diagonal elements of the reduced density matrix go to zero. The Lindblad oper¬ 
ator is 

1 ^ 

^^P) = ^ E - SlSlp - pSlSl) (11) 

k 

where Lfc is the parameter of the kth- qubit which take into account the dephasing of the qubit with the 
environment (the tilde operators do not appear due to commutation of this operators with the evolution 
operator U). 


(d) Correlated dephasing : Here, one takes into account the collective effect of the environment to 
the qubits. Lindblad’s operator is of the form 


^(P) 


1 

ih 


N 


^jk{2S^kPS^ - S^SIp - pS^Sl) 
j,k 


( 12 ) 
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where Tjk = Tkj is the parameter with take into account the correlation ( Ta = Fj). 
The analytical solutions for these four cases are given in the appendix. 


4 Entanglement and GME-concurrence 

Our 3-qubits Hilbert space H is generated by the basis {|CiC 2 C 3 )}?i=o,i- Labeling the qubit of the 3-qubits 
chain as ABC, we understand an entanglement of the form AB when the qubits {1,2} are entangled, and 
we understand an entanglement of the form ABC when the 3-qubits {1,2,3} are entangled. The entangled 
state under our consideration are listen on the table 1. We chose these state since they are mostly used on 
experiments of quantum computation or quantum information. 


Table 1: Entangled states (ordered according their energy separation) 


Entangled form 

Initial entangled state 

AE/{h2TTMHz) 

ABC 

|vhi 8 ) = (|l) + | 8 ))/y 2 ; 
15-27) = (| 2 ) + |7))/y2; 
1^1^36) = (|3) + |6))/y2; 
15/45) = (|4) + |5))/y2; 

= Eg — El = 700 
/S.E 27 = E'j — E 2 = 500 
^EgQ = Eq — Eg = 300 

A _£/45 = Eg — £^4 = 100 

AB 

I«i7) = (|l) + |7))/y2; 

1^28) = (| 2 ) + | 8 ))/y 2 ; 

|a46) = (|4) + |6))/y2; 

l«35) = (|3) + m/y/2- 

AEir = Er- Ei= 605.2 
AE 2 S — Eg — E 2 — 594.8 
AEjig = Eq — £4 = 209.8 
A £35 = £5 - £3 = 195.2 

BC 

|/3i4) = (|l) + |4))/y2; 
1/358) = (|5) + |8))/V2; 
1/323) = (|2) + |3))/V2; 
1/367) = (| 6 ) + |7))/V2; 

A£i 4 = £4 - £1 = 305.2 
A£58 = Eg — Eg = 294.8 
A £23 = £3 — £2 = 104.8 
A £67 = £7 — Eg = 95.2 

AC 

ia 6 ) = (|l) + | 6 ))/V 2 ; 
168) = (|3) + |8))/y2; 
165) = (|2) + |5))/62; 
167) = (|4) + |7))/62; 

A£i 6 = £6 - £1 = 510 
A £38 = Eg — £3 = 490 
A£25 = Eg — £2 = 300 
A £47 = £7 — £4 = 300 
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In order to quantify the entanglement of a state formed by three qubits basis-states, we will use the 
criteria given on [llllliig where the lower bound of the concurrence is 


Cgme{^) > 2 




(13) 


where |$) is a separable state of the two copies of the Hilbert space, Hjo,} is the permutation operator 

acting on the double copies of the Hilbert space, interchanging elements of one space into the other, 

for example 0 IV'iV’ 2 )) = \'ipi4>2) 0 \4>ii’2)- if Ci = {llrn.n}, Ii{iim,n} acting on I'I'rs) means 

that the qubit label by is fixed and the qubits label ”to” and ”n” are interchanged (n{i|2,3}|18) = |45), 
n{i 2,3}|18) = |81) ). Denoting by pmn = (w|'hi8)('hi8|?T-)) one has the fallowing GME-concurrence associated 
to the given entangled state: 


Table 2: Entanglement form: ABC. 

State Id)) GME-concurrence 

I'E'ls) |18) 2 |pi8 | — 2^P44P55 — 2y/p33PeQ — 2^p22P77 

l'E' 27 ) |27) 2|p27| — 2^PiiP88 — “^y/PSSPee — 2-^P44P55 

I'E'ae) |36) 2|p36| — 2^piip88 — ^\/P22Pn — 2-^P44P55 

1 ^ 45 ) |45) 2|p 45| — 2y^iip88 — 2-yp22P77 — 2-yp33P66 


For the case when entangled state is of the form AB, BG or AC, one makes the trace on the missing letter 
qubit, and it follows that ( 1$) = \ijkl) ) 


Cqme 

2 


\{il\pAB\kj)\ - \/{'>'j\PAB\ij){kl)\pAB\kl), 


(14) 


Table 3: Entanglemet form: AB. 


State 

1$) 

GME-concurrence 

0 : 17 ) 

0 : 28 ) 

10101) 

2 Pl7 + P28 

~ (P33 + P44)(P55 + Pee) 

0:46) 

O 35 ) 

10011) 

2 P35 + P46 

~ (Pll + P22)(P77 + Psa) 



Table 4: Entanglement form: BC. 

State 

1$) 

GME-concurrence 

^ 14 ) 

10101) 

2 Pl4 + P58 

— 2-\/ (p22 + P66)(P33 + P 77 ) 

/358) 

/323) 

10011) 

2 P 23 + P67 

~ 2-\/ (pii -|- P55)(P44 + Paa) 

IPer) 
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Table 5: Entanglement form: AC. 


State I $) GME-concurrence 

l^ie) |0101) 2|pi6 + PssI — 2-\/ (p22 + P44)(P55 + Pit) 

l^s) ■ ■ ■ _ 

1 ^ 25 ) |0011) 2|p23 + P 47 I — 2-^/ (pii + /333 )(p 66 + Pss) 

1 ^ 47 ) 


5 Results 


In our case, we have three qubits space {|C 2 '? 2 Ci)}^i=o,i}, and our parameter in units 2t:MHz are 


Wi 

= 400; 

a ;2 = 200 ; 

ws = 100 J = 

10; J' = 0.4 



7i 

= 0.05; 

72 = 0.05; 

73 = 0.05 Ti 

= 0.05; r 2 = 

0.05; 

r 3 = 0.05 

7l2 

= 0.05; 

723 = 0.025 

713 = 0.0125 

Tis = 0.05; 

r23 = 

0.025; ri 3 = 0.0125 


the time is normalized by the same factor of 2ttMHz. To determine the departure of the pure state entangled 
state, we use the purity parameter, P = Tr{p^) [25] . Figure la shows the behavior of this parameter for 
the entangled state I'l'is) as a function of time, where one can see that correlations does not affect much 
the independent model of the environment, which can be seen only for much bigger dissipation parameters. 
Figure lb. Dephasing models finish with the mix state on the system at the end, instead of a pure state of 
the independent model. As seen on Figure 2, independent model ends with a pure state in the system due 
to the system ends on the ground state after sharing energy with the environment. 
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INDEPENDENT 

I'P^^) = (11) + I8))/V2 



Figure 1: Independent Model: (a) Purity behavior with 4-models environments, (b) Indepent model with big 
dissipation. 


I'1'j3) = (I1) + I8)W2 



l'l'j,) = (ll) + l8»W2 


Ti = Ti = t = “5: 1'.2"1'23 = Ti3 = 1» 














concurrence, and Figures 3b, 3d, 4b, and 4d show their associated Purity parameter behavior. The sys¬ 
tem always finish on the pure ground states (|1) = |000)), the lowest bound of the GME-concurrence fall 
down (although this parameter can not tell us whether or not the entanglement has been completely de¬ 
stroyed) . Except for the entangled state Id^is) (maximum entergy difference between their entangled qubits), 
there is not clear difference how this entanglement decay is developed. For example, the entangled states 
I'h 27 ), l^'ae), |'I' 45 ) have the same GME-concurrence decay behavior, but these states have different energy- 
difference on their associated qubits. The entangled states \an) and |Qf 46 ) the GME-concurrence decay is 
the same, although their energy-difference is quite big. 




(a) (b) 

Figure 3: Independent:(a) GME-concurrence vs t. (b) Purity vs t. 
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Figure 4: Independent:(a) GME-concurrence vs t. (b) Purity vs t. 
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Dephasing Model: In this case, the GME-concurrence and purity parameters can be expressed explicitly 
in a simple form as shown in the next table: 


Table 6: 


Entanglement 

CGME-.\rii){T) 

P\r,,){r) 

{ABC) 

{AB) 

{BC) 

{AC) 

2|p„(0)exp(-rT)|, 
2|p„(0)exp[-(ri -br 2 )T]|, 
2|py(o) exp[-(r2 -br3)T]|, 
2|p„(o)exp[-{ri -br3)T]|, 

PiM + pjj(o) + 2py (0) exp(-rT) 
pfi(O) + pi(0) + 2p|(0) exp[-(ri -b r 2 )T] 

pfi(o) + pi(o) -b 2 p|(o) exp[-(r 2 -b r3)T] 

pfi(O) + P^j(O) + 2pfj-(0) exp[-(ri -b r 3 )T] 


where T = Ti +r 2 + r 3 , and |Ty) = {|d>ij), \aij), |/3y), for the entangled cases ABC, AB, BC and AC 
respectively. This expressions show that the decay behavior is the same for each family of entangled states, 
that is, entangled states in the same family have the same decay behavior. 


Independent correlated: From Figure 1, we saw that correlations have not effect on the purity. In 
addition, Figure 5 shows entangled states in different environments where we see that the behavior of the 
GME-concurrence is the same for the independent and independent correlated models. 



(i)3 



(i)3 


(a) 


(b) 


Figure 5: GME-concurrence for entangled states with the four difference environments. 
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6 Conclusions 


We have studied the decay behavior of entangled states, formed by two basic states of three qubits registers, 
under four difference environments and using Lindblad type of equation to see whether or not this decay has 
a regular dependence with respect the energy-difference (difference of energy of the two basic states of three 
qubits which made up the entangled state) associated to the entangled state. We did not find this regular 
dependence, but rather a complicated situation which depends also on the type of environment. 
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Appendix 


We consider a linear chain of three nuclear spin system. Then, our basis is {|? 3 C 2 ^i)} 5 i=o,i) and the equations 
for the reduced matrix elements are obtained by making the bracket with these elements of the basis of the 
equation 


(a) Independent : 

-^Pmn + ^ ^5mn{k)pmn = p^+2^-k ^^+2^ ( 15 ) 

k k 

where we have made the following definitions 



= 5aT- i5 m- 2 N-fc +l5a",l(5 „_2N-fc n 

(16) 


= ^a^,0^a^,0 

(17) 


^fc,n+2-^“^ 771+2-^“^ ■ 

(18) 

(b) Correlated independent : 




_5 

Pmn\t) 


Y, [26muik, 


k,l 


-Smik, 

-dn(K,l)e Pm,n-2^-^+2^-fc 


(19) 


where the following definitions have been made 


— SaV^^lS Tn-2^-^ n 

i ,u 



(20) 

- ^l,m-2^-‘ 

(21) 

- llz,n-2W-l 

(22) 


(c) Dephasing : 


N 


dt 


Pmn(t)=^rfc (-1) 


+“fc _ X 


Pmn (t) 


(23) 
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which has the following analytical solution 


Pmnii^ — Pmn(O) OXp 


N 


i-(-i) 




(24) 


(d) Correlated depahsing : 


df 


.it) = 


N 

E 

k,l 


^kl 

~T 


2(-l) 




-(- 1 ) 




-(- 1 ) 






which has the explicit solution 


(25) 


Pmn(i) — Pmn(O) OXp \ 



(- 1 ) 






2(-l)^ 


(26) 
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